OUT 4102 54 



hep-ph/9504398 



INP 95 13/377 



Three— loop QED Vacuum Polarization and the 
Four— loop Muon Anomalous Magnetic Moment*^ 



p. A. BAIKOV* 

Institute of Nuclear Physics, Moscow State University 
119 899, Moscow, Russia 

and 

D. J. BROADHURST§ 
Physics Department, Open University 
Milton Keynes, MK7 6AA, UK 



Three-loop contributions to massive QED vacuum polarization are evaluated by a com- 
bination of analytical and numerical techniques. The first three Taylor coefficients, at 
small q-^, are obtained analytically, using d— dimensional recurrence relations. Combining 
these with analytical input at threshold, and at large q^, an accurate Pade approximation 
is obtained, for all q^. Inserting this in the one-loop diagram for the muon anomalous 
magnetic moment, we find reasonable agreement with four-loop, single— electron-loop, 
muon-anomaly contributions, recently re— evaluated by Kinoshita, using 8— dimensional 
Monte— Carlo integration. We believe that our new method is at least two orders of 
magnitude more accurate than the Monte-Carlo approach, whose uncertainties appear 
to have been underestimated, by a factor of 6. 

1. Introduction 

We describe a metltod, previously tested^ in two-loop QCD, to approximate, to high 
accuracy, three-loop contributions to QED vacuum polarization, using new an; 
ical results foLthe small momentum-transfer limit, combined with asymptotici 
and thresholdO results. Related contributions to the four-loop muon anomalous 
magnetic moment are computed, to test an evaluationEl that was undertaken in 
response to a previous discrepancy between numericalQ and analyticaltl work. 

In the on-shell (OS) renormalization scheme of conventional QED, the renor- 
malized photon propagator has a denominator (1 + n(z)), where z = /Aw?, 
with an electron mass m, and the vacuum polarization function, n(z), vanishes at 
z — Q. Non-relativistic consideration of the electron-positron system yields infor- 
mation about the thresholdd behaviour, as 2 ^ 1. Moreover, the MS asymptotic 
behaviourSa, as z —00, combined with relationsBa between the MS and OS 
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schemes, yields two terms of the asymptotic expansion in powers of 1/z. 

The crucial new ingredient is our use of recurrence relationsLi, to obtain the 
three terms of the expansion as z — > 0. Combining 6 analytical data with Padc 
(or hypergeometridd'EJ) approximations, we shall produce reliable fits, for all z, and 
hence check four-loop muon-anomaly contributiona3. 

2. Small momentum behaviour 

We evaluated, to 3 loops, the first 3 moments in the z ^ expansion 

n(z) = ^C„z" + 0(a4), 

n>0 

by intensive application of d-dimensional recurrence relations to three-loop mas- 
sive vacuum diagramsa, with propagators raised to powers up to 11, since up to 
8 differentiations w.r.t. the external momentum q are required before^ setting it to 
zero. This put great demands on the REDUCE package RECURSORB, which used 
80 MB of memory, for 2 days, on a DecAlpha machine, after hand-tuning the pro- 
cedures, to minimize recomputation of integrals, and to allow safe truncation in 
e — (A^~d)/2. The gauge invariance of Ci and C2 was verified for all e. After OS 
massLJilj and chargeQ renormalization, we obtained the finite e ^ limits 
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where wc follow common practiced, by allowing for N degenerate leptons. In pure 
QED, = 1; formally, the powers of N serve to count the number of electron loops. 
Our principal interest, for consideration of four-loop muon-anomaly contributions^, 
is Hg^'a^/TT^, the three-loop contributions to 11 that involve a single electron loop. 
The moments of n3^'(z) are given by the coefficients of Na^/ir^. (We shall not 
need the N'^a^ /ir^ terms in Section 6, since the muon-annmaly contributions of the 
diagrams with two electron loops are better understoodBEJ.) 
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3. Large momentum behaviour 

The situation regarding the z —> —oo behaiiiour of n(z) was unclear, until re- 
cently, becajise three-loop MS QCD resultsB had been altered, while obtaining 
QED resultaJ, in the belief [now knownP to be mistaken) that the former contained 
errors. Further calculationEJ confirmed the QCD resultscl and hence invalidated 
the 0(l/z) QED resultsB. Accordingly, we thought it prudent to derive the OS 
asymptotic behaviour ourselves, from first principles, using the REDUCE package 
SLICER, which hacLbeen written specifically to checkB the leading, massless, MS 
behaviour, obtained! with the SCHOONSCHIP package MINCER0. 

In our ab initio derivation of the asymptotic OS result for Ilg^' , we used neither 
the MS scheme, nor MINCER. Instead, the asymptotic expansion of the bare dia- 
grams was obtained, in d dimensions, using SLICER, and the bare charge and mass 
were transformed. liirectly to the physical charge and mass, using multiplicative OS 
renormalizationsO'Ej, obtained by RECURSOR. Setting e — 0, we obtained a finite 
OS result of the form 

nW(z) = A{z) + B{z)/z + 0{L^/z^) , 
where L = ln(— 4z) = \n{—q^/m?) and 

= -Ti + iC5-iC3 + 2(ln2-|)C2 + ^L, (4) 
B{z) = ^-iC5-iC3 + 3(ln2-|)C2-i(i-6L2) . (5) 

Using finite transformationsBlli from physical to MS-renormalized quantities, one 
obtains, from our OS result, an MS asymptotic behaviour identical to that which 
wouldEZi have been obtained from the QCD analysisB, had it not been niiscorrected 
in the course of deriving QED resultsQ. As a result of our, and otheilij, work, a 
(second) erraturm to the QED workS was issued. 

In conclusion, we are con fident of our OSJDED result, since it is quite indepen- 
dent of previous workaj'BEJ and, eventuallyuO, in agreement with them. 

4. Threshold behaviour 

The leading threshold behaviour, at 3 loops, is determined by non-relativistic quan- 
tum mechanics: n3^'(2;) = ^""^(l ~ z)^^/'^ + 0(ln(l — z)), as z ^ 1. Moreover, it 
appearsu'113 that a stronger statement can be made, namely that the first relativistic 
correction to the spectral density, p{t) = Imn(t -I- i0)/7r, at any given order in a, 
is cancelled in the combination (1 -I- Aa/n)p{t). At the two-loop level, the exact 
relativistic results confirm that P2{t) + 4 pi{t) = tt^ + O(u^) is free of a term of firet 
order in u = (1 — 1/t)^/^. The corresponding cancellation at 3 loops is expectedlla 
to occur in v{p^^^ + 4p2) = ^i''* + 0(^^^); implying that 

lim fn[,^'(z) + 4n2(z) ^ \ = constant , (6) 



24(l-z)i/2 



with an unknown value for the constant, but no logarithmic singularity. 
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5. Approximation method 

We express the analytical results as properties of the combination 

nw (z) ^ nw (z) + 4 (z) + (1 - .) g(z) f ? g(z) + + 



16 32z J 32z 96 

where G(2:) = 2F1 (l, 1; 1;^:) is given by (z^ — z)^^/'^arcsinh(— z)^/^, on the nega- 
tive real axis, and the two-loop term, 112 (z), is quadratic in G(z) and involves a 
trilogarithmB and its derivative. 

The data are conveniently encoded by the moments 

of the spectral density of ftg^'. At small z, we have il\l\z) — X]r!>o^^("') -^^ ^^'^ 
hence obtain M{1), M(2), M{3) from the coefficients of Na^/n^ in the results of 
Eqs (0J^,^) for the moments of n(z), after taking account of the known moments of 
the additional terms in Eq (^. At large z, the logarithmic singularities of these ad- 
ditional terms cancel, by deliberate construction, those of Eqs (^,^ , whose constant 
terms therefore determine M(0) and M(— 1), respectively. Finally, the threshold 
Coulomb singularity of Eq gives the large-n behaviour of the ratio 

RH - ^ = ^ + 0(l/n) , G(n) . I ^ - {rj^ , 

where G(n) is the moment of a spectral density p^^\t) = t~^^'^{t — 1)^^/^, with a 
coulombic 1/v threshold singularity and the same convergence properties, at large 
as Pg^' {t). Note that a further datum, namely the absence of a logarithmic singularity 
in Eq (^), corresponds to the absence of an 0(l/rt^/^) term in R{n), as n ^ 00, 
partly accounting for the remarkable uniformity of our final analytical database: 

R{-1)= -i(ln2-|)C2 + iC3 + iC5 + ii = 3.473721028889 

i?(0)= -i(ln2-|)C2 + ^C3-fC5 + ii = 4.087 577 319 347 

i?(l)= -|(ln2-|)C2 + ™C3 + fii = 4.188 975 919 282 

R{2)= -|(ln2-f)C2 + «fC3-^^ = 4.216 954 083 059 

i?(3)= -I (ln2 - I) C2 + mmWCs - %ttlZlfol' = 4.224 481581719 

i?((X)) = = 4.058 712126 417 

which has been obtained from 3 quite disparate regimes. 

Now we map0 the cut z-plane onto the unit disk, and define a mapped function 

^<"'-(^ '-jrfijr- <«) 

which is analytic for |w| < 1, with the cut mapped to the unit circle. The 6 data 
then determine {P(-l), P(0), P'(OL P"(0), P"'(0), P(l)}, allowing us to construct 
[2/3] and [3/2] Fade approximantsEJ, with benign poles outside the unit disk, and 
imaginary parts on the unit circle that accurately approximate the spectral density. 
The differences between these two approximations are very small, for all < 1. 
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6. Four loop contribution to the Muon Anomalous Magnetic Moment 

Our simple rational approximations to P(lu) reproduce, exactly, all known data on 
n]^^ , as well as its analyticity structure. We now use them to calculate the four-loop 
contributionl§, = A^^'a^/Tr^, to the muon anomaly, {g/2 — 1)^, due to insertion 
of three-loop, single-electron-loop vacuum polarization diagrams into the one-loop 
anomaly diagram. A typical diagram is 




The resulting coefficient of a'^/n^ in the muon anomaly is given b; 



4' = 



dy(l-y)nW(z) 



4™^ I — y 



(9) 



We calculate the integral using [3/2], [2/3], and [2/2] Pade approximants to P{lo)- 
In the [2/2] approximants we omit a piece of data from each regime, obtaining 



Input 
Pade 



-A, 



[1] 



all 
[3/2] 
0.230 362 20 



all 
[2/3] 
0.230 36218 



omit Eq (|) 

[2,2] 
0.230 360 42 



omit Eq 
[2/2] 
0.230 366 94 



omit Eq 
[2/2] 
0.230 36149 



with a muon mass = 206.768 262 m. The stability is remarkable: changing the 
Pade method from [3/2] to [2/3] changes the output by 1 part in lO''; removing 
a piece of data, from any of the 3 regimes, changes it by no more than 2 parts 
in 10^. The improvement from using 6 inputs, as opposed to 5, is greatest in the 
case of including the asymptotic result of Eq (H). In contrast, the Coulomb datum. 



R{oo) 



improves the result by only 3 parts in 10 , since the muon-anomaly 



integral (Q) involves only space-like momenta. The smallness of our spread of results 
demonstrates a high degree of consistency in the input, making the possibility of 
analytical error very remote. Being conservative, we take the range of [2/2] results 
as a measure of our uncertainty, and arrive at = —0.230362(5), to be compared 
with a recen to Monte-Carlo result, j^]} = -0.2415(19), obtained using VECAS, in 
preference to RIWIAD (which gaveQ a grossly discrepant value, amendedQ in the 
light of a renormalization-group analysiaS). In visual terms, the comparison is 
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To verify that this (discrepancy is not an artifact of the Pade method, we also 
tried a hypergeomctricIO method, i.e. a polynomial fit to p\i\t) / p^^\t) — ^tt** = 
Efc>i Cfe w^ with 5 coefficients, fixed by R{n) - ^n* = J^k Cfe {1)^+2 I (^)„+2- 
As might be expected, the resultant fit to the spectral density (t) was less smooth 
than in our Pade methods. Nevertheless, the shift in the value for A^^^ was two 
orders of magnitude less than the disagreement with the Monte-Carlo result. 

In conclusion, we stress that the analytical data of Eqs (1-6) exhibit a high 
degree of internal consistency, making it most unlikely that any of them is in error. 
Pade approximants for the mapping (||) of the well-behaved function (|^) enable us 
to evaluate the muon-anomaly contribution with an uncertainty of 2 parts in 
10^. Our result is in reasonable agreement with a recent, lower-precision, Monte- 
Carlo re-evaluationO, whose uncertainties appear to have been underestimated by 
a factor of 6, which is a great improvement on fhe situation revealed by a previous 
discrepancy between analyticalcl and numericaB work. 
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